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Abstract 



Using spectral decomposition techniques and singular perturbation theory, we develop a 
systematic method to approximate the prices of a variety of European and path-dependent 
options in a fast mean-reverting stochastic volatility setting. Our method is shown to be 
equivalent to those developed in [10], but has the advantage of being able to price options for 
P^ ■ which the methods of [10] are unsuitable. In particular, we are able to price double-barrier 

f r^ I options. To our knowledge, this is the first time that double-barrier options have been priced 

in a stochastic volatility setting in which the Brownian motions driving the stock and volatility 
are correlated. 



> ; 1 Introduction 

^-H ! 

^O ■ Since it was originally analyzed in the context of Sturm-Liouville operators, spectral theory has 

^\^ [ enjoyed wide popularity in both science and engineering. In physics, for example, the stationary- 

^. ■ state wave functions are simply the eigenfunctions of the time- independent Schrodinger equation. 

C^ I And, electrical engineers are well-versed in the theory of Fourier series and Fourier transforms. It 

is not surprising then, that techniques from spectral theory have found their place in finance as 
well. 

For instance, in [23] eigenfunction methods are used to price European-style options in a 
Black-Scholes setting. The authors of [16] use eigenfunction techniques in the context of bond 
^^ ■ pricing. Spectral decomposition techniques have been particularly successful at aiding in the 

?H I development of analytic pricing formulas for a variety of exotic options. For example, in [30], 

Fourier series methods are used to obtain closed-form expressions for prices of double-barrier 
options in the Black-Scholes setting. And in [8,25,26] spectral decomposition techniques are used 
to obtain analytic option prices-both European and path-dependent-where the underlying and 
short rate are controlled by a one-dimensional diffusion. Additionally, the authors of [27] use 
spectral methods to evaluate both bonds and options in a unified credit-equity framework. 

Like spectral theory, stochastic volatility models have become an indispensable tool in math- 
ematical finance. By and large, this is due to the fact that two of the earliest and most well- 
known stochastic volatility models-the Heston model [19] and Hull- White model [22]-capture the 
most salient features of the implied volatility surface while at the same time maintaining analytic 
tractability. Stochastic volatility models have become so popular, in fact, that entire books have 
been written on the subject [10,15,24]. 

It seems natural, then, to try to employ elements from spectral theory in a stochastic volatility 
setting. Yet, other than the spectral decomposition of various volatility processes, which is expertly 
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done in [24], there is a surprising lack of literature in this area. In particular, we are unaware of 
any literature that uses spectral methods to price double-barrier options in a stochastic volatility 
setting in which the Brownian motions driving the stock and volatility are correlated. The difficulty 
with using spectral analysis when the stock price and volatility are correlated arises because 
spectral decomposition techniques work best when there is some sort of symmetry inherent in the 
problem being studied. This symmetry is broken when the stock price and volatility are correlated 
via two Brownian motions. Yet we know that correlation between the stock price and the volatility 
is important because it is needed in order to capture the skew of the implied volatility at the money 
and reflect the leverage effect [10, 15]. 

In this paper, we apply techniques from spectral theory to a class of fast mean-reverting 
stochastic volatility models in which the stock price and volatility are correlated via two Brownian 
motions. The two-dimensional diffusion that controls the stock and volatility is in contrast to the 
work of [7,8,25,26,28], where spectral and probabilistic methods are used to price options on scalar 
diffusions. Extensions to two-dimensions are highly non-trivial and it is this that distinguishes 
our work from the earlier contributions. The class of fast mean-reverting stochastic volatility 
models, first studied in [10], is an important class of models to consider because volatility has 
been empirically shown to operate on short time-scales [11,20]. To price options in this setting 
we employ the singular perturbation methods of [10], but we do this in the context of a spectral 
expansion. 

The rest of this paper proceeds as follows. In section 2.1 we introduce the class of fast mean- 
reverting stochastic volatility models first considered in [10]. Additionally, we discuss how this 
class this class of models relates to two of the more popular models used in practice - SABR and 
Heston. In section 2.2 we present an option-pricing framework, which allows us to consider both 
European, single- and double-barrier options. This framework results in an option-pricing partial 
differential equation (PDE) along with appropriate boundary conditions (BC's), which must be 
solved in order to specify the price of an option. We briefly mention how the authors of [10] use 
singular perturbation theory to obtain an approximate solution to the option-pricing PDE and 
explain why for certain options (e.g. double-barrier options) this methodology is unsuitable. We 
then present a new method of solving the option-pricing PDE - one which is suitable in cases 
where the methods of [10] are not. In this new method we assume a solution of a speciflc form and 
show how this leads to an eigenvalue equation. An approximate solution to the eigenvalue equation 
is given in section 3.1. Then, using this solution, we provide formulas for the approximate price 
of an option in section 3.2. Equivalence of the option-pricing formulas presented in this paper to 
those derived in [10] is established in section 3.3, as is the accuracy of our pricing approximation. 
In section 4 we discuss the practical implementation of our methods. We present three examples: 
European calls, up-and-out calls, and double-barrier knock-out calls. Additionally, we mention 
how our framework can be extended to price knock-in and rebate options. We flnish by discussing 
some issues related to calibration. 

2 A Class of Models and an Outline of Our Method 

In this section we introduce a class of fast mean-reverting stochastic volatility models. We then 
present an option-pricing problem and outline our method for obtaining an approximate solution 
to this problem. 

2.1 A Class of Fast Mean- Reverting Stochastic Volatility Models 

We study the class of fast mean-reverting stochastic volatility models flrst considered by Fouque 
et al. in [10]. Specifically, under the risk-neutral pricing measure P, we consider a non-dividend 
paying asset (stock, index, etc.) St = exp (Xt) whose dynamics are given by the following system 



of stochastic differential equations (SDE's 

1 



dXt^[fi--f{Y,^)]dt + f{YndWt 
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d{W,B)^ ^ pdt. 

Here, Wt and Bt are Brownian motions under P with instantaneous correlation p such that p^ < 1. 
The price process St follows a geometric Brownian motion with growth rate /x, which equals to 
the risk- free rate of interest, and with stochastic volatility fiY^) (the traditional symbol for the 
risk- free rate of interest r is reserved for a different purpose). The dynamics of Xt — \og St are 
obtained from Ito's Lemma. We note that, as it should be, the discounted stock price (e"^*S't) 
is a martingale with respect to the canonical filtration of the Brownian motions. The process Y^ 
evolves as an Ornstein-Uhlenbeck (OU) process under the physical measure P. That is 

dY^ — - (y -- Yf) dt ~i r=r- dBt. (under physical measure P) 

e Ve 

However, under the risk-neutral measure the dynamics of Y^ acquire a market price of volatility 
risk, which is given by K{Y^). The superscript on Y^ indicates that this process evolves on a time- 
scale e. The parameter e is intended to be small (i.e. < e <C 1) so that the rate of mean-reversion 
of the OU process is large. It is in this sense that the volatility is fast mean-reverting. We note 
that under the physical measure, Y^ has a unique invariant distribution Y^ ~ A/'(y, u^). 

It is not necessary to specify the precise form of f{y) or A(j/), as only certain moment will 
play a role in the asymptotic analysis that follows. Likewise, the particular choice of Y^ as an OU 
process is not crucial for our analysis. However, in order to guarantee the accuracy of our pricing 
approximation we need the following assumptions: 

1. Under the physical measure, Y^ has a unique invariant distribution, which is independent of 
e. 

2. Under the physical measure, the moments of Y^ are uniformly bounded in t. Note that this 
assumption actually follows from the previous assumption. 

3. The smallest non-zero eigenvalue of Cy - the infinitesimal generator of Y^ under the physical 
measure - is strictly positive. 

4. There exists a constant Ca > such that |A(2/)| < Ca. 

5. The function /(y) is such that the solution (t){y) of Poisson equation (3.5) is at most poly- 
nomially growing. 

We note that the Cox-IngersoU-Ross (CIR) process, as well as the OU process satisfy the above 
assumptions. 

Of practical interest is how two of the most popular stochastic volatility models - Heston and 
SABR - fit within the fast mean-reverting class of models discussed in this paper. The Heston 
model [19] can be accommodated in this class by choosing Y^ to be a CIR process under the 
physical measure, choosing f{y) — ^/y, and setting A(y) — 0. The rate of mean-reversion of the 
CIR process should then be scaled by 1/e and the "volatility of volatility" term should be scaled 
by l/-\/e. This will ensure that the invariant distribution is independent of e. Note that the choice 
A(y) := is not really a restriction on the Heston model as the stochastic variance Y^ in Heston 
is a CIR process under both the physical and risk-neutral measures. 

The key change between the Black-Scholcs model [3] and the Heston model is that the constant 
volatility of Black-Scholes is replaced by the square-root of a CIR process in Heston. In this sense. 



the Heston model is essentially an extension of the Black-Scholes model. As such, the CIR process 
that controls the volatility in the Heston model can be treated as a perturbation around geometric 
Brownian motion with constant volatility, just as the fast mean- reverting factor of volatility is 
treated as a perturbation around constant volatility geometric Brownian motion in this paper. 
This is the reason that the Heston model fits within the fast mean-reverting stochastic volatility 
framework. 

Contrary to Heston, the SABR model [17] is an extension of the Constant-Elastic- Variance 
(CEV) model [6] . As such, the SABR model does not fit within the class of models considered in 
this paper. That said, an extension of the CEV model that includes a fast mean-reverting factor 
of volatility is possible. Approximate option prices for the CEV model with a fast mean-reverting 
factor of volatility are derived using singular perturbation theory in [5]. Also note, because option 
prices in the CEV framework have a spectral representation [8] , the combined singular perturbation 
and spectral method developed in this paper would be suitable for CEV with a fast mean-reverting 
factor of volatility. However, this is outside the scope of the present work. 



2.2 Statement of the Option Pricing Problem and Outline of our Method 

In this section we introduce an option-pricing problem and outline our method of obtaining an 
approximate solution to this problem. 

Consider an option expiring at time T < oo whose payoff can be expressed 



h: lU{l}U{r} 



T = inf{t > : Vt ^ /} A T, 
— oo < I < r < oo, 
h{l) = h{r) = 0, 



(2.1) 



In words, the option has payoff h^Xx) if X does not leave the interval / prior to time T, otherwise 
the option payoff is zero. We use the convention inf {0} — oo. Options that fit within the above 
payoff framework include European and knock-out style options. But, we shall see in section 4.4 
that our results can be extended to include knock-in and rebate style options as well. 

We introduce a money market account Mt = e^*, which we shall use as our option-pricing 
numeraire. According to risk-neutral pricing, the value of the option P^ at time s < T is given by 






h{Xr) 



Mr 



Ts 



J-, = a({Xt,y/:0<t<s}), 



After a bit of algebra, and using the Markov property of [X^V^] one finds 



Pt = \r<s}e^^'-^^h{Xr) + E [e-^^(--'^h{Xr) \r>s} 
= I{.<.}e^(^-^)/i (Xr) + l{r>s}P' (s, V„ y/) , 



Ts 



where 



and 



F^(s,X„r/)=Efe-^(-°-^)/i(X,e) 



x.,y: 



(2.2) 
(2.3) 



r" = inf{t >s:Xt(^I}AT. 

We note that the first term in (2.2) is zero as X^ is either equal to I or r on the set {r < s} and 
h{l) = h(r) = 0. Thus, the price of the option at time s is given simply by the second term of 
(2.2). From (2.3), one deduces that the function P'^{s,x,y) solves the following PDE and BC's 
(see e.g. Chapter 9 of [29]) 

0^ {ds-fi + £x.Y)P', {s,x,y)e [0,T]x/xR, 
h{x)^P^T,x,y), 

= P"(s,?,2/), ff />-oo, (2.4) 

= P"(s,r,2/), ff r < oo. (2.5) 



Note that BC's (2.4) and (2.5) are not required if Z or r are infinite. The notation L\ y represents 
the infinitesimal generator of (X, y^). For clarity, we write CfxY explicitly and state its domain 
V 

^x.Y = \{{y- y) dy + v'dl) + -i= [pvV2f{y)dly - vV2A{y)dy) 

-l.fiy))d. + lfiy)dL, 

V = I xR. 

To simplify subsequent calculations we introduce u''-{t,x,y) such that 

P'is,x,y) = e-'''u'it,x,y), t = T-s. (2.6) 

A straightforward substitution shows that u^{t,x,y) satisfies the following PDE and BC's 

0={-dt+£'x^Y)u", it,x,y)e[0,T]xIxR, (2.7) 

hix)^u'iO,x,y), (2.8) 

= u%t,l,y) if l>-oo, (2.9) 

0^u''{t,r,y) if r < oo. (2.10) 

We set u'^{t, x,y) = for x ^ I. Although u'^{t, x, y) is in fact the un-discounted price of an option 
with time-to-maturity t — T — s, from this point onward we shall refer to u'^{t,x,y) simply as 
the price. For convenience, the theorems derived in section 3.2 are given in terms of u'^{t, x, y), as 
are the examples provided in section 4. The reader should keep in mind that the true price of an 
option P^{s,x,y) at time s can be recovered from u'^(t,x,y) using (2.6). 

In [10] the authors use singular perturbation techniques to find an approximate solution to 
PDE (2.7) by expanding u'^{t, x, y) in powers of the small parameter ^/e 

Roughly speaking, the authors of [10] show 

1. The functions u^^> and u^^> are independent of y. 

2. The the O (e") price is given by u^^\t^x) — u^^{t,x), where u^^{t,x) is the Black-Scholes 
price of an option (with an appropriate level of volatility) . 

3. The the O (e^/^) price u^^) solves C^'^u^^'> = A^^'>u^'^, where C^'^ ^ {-dt+Cx) is the 
Black-Scholes pricing operator, Cx is defined in (3.32) and A^^^ is a linear operator defined 
in (3.6). 

For European-style options - for which u'^{t,x,y) must satisfy only BC (2.8) - and for single- 
barrier options - for which u'^(t,x,y) must satisfy only BC's (2.8) and one of either (2.9) or 
(2.10) ~ the method of [10] works well because in these cases there exist analytic formulas for 
u^ {t,x). However, for double-barrier options - for which u'^{t,x,y) must satisfy all three BC's 
(2.8), (2.9) and (2.10) — the methods of [10] are problematic because the Black-Scholes price of 
a double-barrier option must be expressed as an infinite series [30]. 

In this paper, we use a combination of singular perturbation techniques and spectral methods 
to solve PDE (2.7) with BC's (2.8) - (2.10). The spectral method is outlined as follows: suppose 
we have the solution to the following eigenvalue equation 

J^x,Y^l = KK (2-11) 

= *^(?,2/) if l>-^, (2.12) 

= *^(r,y) if r<oo. (2.13) 



By "solution to the eigenvalue equation" we mean that we have the full set of eigenvalues A^ and 
corresponding eigenfunctions '^Ux, y) for which which equations (2.11) - (2.13) are satisfied. Then 
it is clear that any linear combination of functions of the form e''*9*^^(x, y) will satisfy PDE (2.7) 
and BC's (2.9) and (2.10). Hence, as long as the eigenfunctions allow us enough flexibility to 
match EC (2.8), the function u'^{t,x,y) can be expressed as ^ 

u^-it, X, y) = [^fij^l^f'^^'l f^'T^ ^^^^'^7^ , , m . exp (A^t) ,(2.14) 

\ j A]^ gl[t) ^]^[x T y) dv (contmuous spectrum) ^ -■ h ' 

where q is a place-holder for either n oi v and A'L are constants to be determined by the payoff. 

The main advantage of the spectral method is that by separating the spatial variables (x, y) 
from the temporal variable t the BC's (2.9) and (2.10) can be dealt with at the level of the 
eigenfunctions \E'^(x,j/) as in (2.12) and (2.13), rather than at the level of option prices u'^{t,x,y) 
as in (2.9) and (2.10). This method is particularly advantageous for pricing double-barrier options. 

Using representation (2.14), what remains in order to specify the price of an option is to find 
expressions for \l/^(a:,y), A^ and A^. This is the subject of section 3. 

3 Asymptotic Analysis 

An outline of the asymptotic analysis performed in this section is as follows. First, in section 3.1 we 
derive an approximate solution to eigenvalue equation (2.11). The key results of this derivation are 
presented in Theorem 3.1, Proposition 3.2 and Theorem 3.3. Next, in section 3.2 we use the results 
of section 3.1 to derive an expression for the approximate price of an option. This expression, which 
is given explicitly in Theorem 3.4, serves as the main result of this paper. In section 3.3, we prove 
that our method of obtaining the approximate price of an option is equivalent to the method 
of [10]. We summarize this equivalence in Theorems 3.6 and 3.7. Finally, in Theorem 3.8 we also 
establish the accuracy of our pricing approximation for the case of European options. 

3.1 Asymptotic Solution to the Eigenvalue Equation Cxy^l = K^g 

For general /(y) and A{y) there is no analytic solution to the eigenvalue equation £%^ y'^q — ^q'^q- 
However, we note that C)^ y can be conveniently decomposed in powers of ^/e as 

C^-'^ = iy-y)dy + v'dly, 

£(-1) = pvV2f{y)dly - vV2A{y)dy, 

>C(")-(/i-^f(y))9. + ^/^(y)aL. 

This decomposition suggests a singular perturbative approach. To this end, we expand ^^ and A^ 
in powers of ^/e. We have 

A^ = A(0) + ViA«+6Af) + .... 



^For simplicity wc have assumed either a purely discrete or absolutely continuous spectrum. In fact, depending 
on the operator £^ y and the BC's, the spectrum may be discrete, absolutely continuous or mixed. However, in this 
paper we do not endeavor to solve the full eigenvalue problem (2.11)-(2.13). Rather, we use singular perturbation 
techniques to find an approximate solution to (2. 11)- (2. 13). For the asymptotic analysis we perform in section 3.1 
we shall need to consider only discrete or continuous spectra. 



We now insert the expansions for ^^(a;, y) and A^ into eigenvalue equation (2.11) and collect terms 

of like-powers of -^e. The 0{e^^) and 0{e^^^^) equations are 

O(e-i): = £(-2)^(0)^ 

0(e-i/2). = £(-2)^^(1) +^(-1)^(0)^ 

Noting that all terms in C'^^^'' and C^^^'' take derivatives with respect to y, we may choose 
solutions of the form ^q — ^q (x) and ^q — ^^ (x) (i.e. functions of a; only). Continuing the 
asymptotic analysis, the order ©(e") and ©(e^/^) equations are 

O(e0) : £(-2)^(2) ^ (^^(o) _ ^(o)^ ^(o)^ (3^^) 

0(ei/2) : £(-2)vI/(3) = -£(-i)M/(2) + (a^^' - £W) *« + A«vl/(o), (3.2) 

where we have used C'^^^^'^q (x) — in (3.1). Equations (3.1) and (3.2) are Poisson equations for 
^q {x,y) and ^q {x,y) respectively in the variable y with respect to the operator £(^^' = eCy. 
We remind the reader that the operator e Cy is the infinitesimal generator of Y^^ under the physical 
measure. In order for an equation of the form CYu(y) ~ v{y) to have a solution with reasonable 
growth at infinity, the following centering condition must hold 

v{y)dFY{y)^Q, 

where Fy is the invariant distribution of Y^ under the physical measure. For the prototype OU 
used in this paper Fy ^ A/'(y, w^). Throughout this paper, the notation (•) will always indicate 
averaging with respect to the invariant distribution Fy . In equations (3.1) and (3.2) the centering 
conditions become 

0=(a(°)-(£(°)))*(°), (3.3) 

- - (£(-i)*f ) + (a(o) - (£(0))) vI/W + A«*(o). (3.4) 

Using appropriate BC's, eigenvalue equation (3.3) can be solved exphcitly, as the operator (£('^') 
is given by 

However, in order to solve eigenvalue (3.4), we need an expression for (£(^^)^^ (a;, •) V To this 
end, we note from (3.1) 

£(-2)*f) ^ (a(°) - £(°)) *(°) ^ ((£(")) - £(")) *(o) ^l{.'- f) (9L - a.) *("). 

Now, introducing 4>{y) as a solution to the following Poisson equation 

£(-2)0 = /2-a2, (3.5) 

wc may express Wq {x^ y) as 

'^^^\^.v) = -\c^{y){dL-d.)^^:;\x). 
Hence, (C^^^^^q (x, ■)) is given by 

L^-'^-^f) = UpvV2f{y)dly vV2k{y)dy) (-\<j,{y) (9^ - 9.) ¥^Hx)) ) = ^^^^ *f (x). 



where 

A^'^ = V, (5L. - aL) + V, (aL - 9.) , ^2 - ^ (a^') , ^3 = ^ {l<i>') ■ (3.6) 

Thus, from (3.4) we have 

^(1) vi/(o) ^ (aW - (£(0))) *« + A«*(o). (3.7) 

Given a solution to (3.3), one can use (3.7) to find expressions for Wq ' {x) and \q . 

FinaUy, we make a remark about BC's. In order to satisfy (2.12) and (2.13), we must also 
impose the following BC's 

$(o)(/) = ^(i)(/)^0 if />-oo, (3.8) 

¥^\r) = *^^'(^) =0 if r < oo. (3.9) 

This is as far as we shall take the asymptotic analysis. The key results of this analysis equations 
(3.3) and (3.7) , which can be used to find expressions for ^\ {x), Xq and ^g '{x),Xq ' . We shall 
present these expressions in Proposition 3.2 and Theorem 3.3. Before doing so, however, we 
establish some key facts about the eigenfunctions ^^ '{x) of (3.3). 

Theorem 3.1. The eigenfunctions ^q (x) of equation (3.3) form a complete orthonormal basis 
in the Hilhert space % := L^{I, s) where 

s{x) dx — —e '^^dx, c= ^ , {u,v)^— / u{x)v{x)s{x)dx. 

o" '^ Ji 

In the case where both I and r are finite, the spectrum is discrete. In all other cases the spectrum 
is absolutely continuous with respect to the Lebesgue measure. 

Theorem 3.1 is a standard result from Sturm-Liouville theory and can be found in any number 
of texts on differential equations [1,2,21,31,33]. A sketch of the proof of Theorem 3.1 is as follows. 
First, note that (3.3) may be recast in standard Sturm-Liouville form 

9x (e^^^a,*^ (x)) = AW six) ¥^Hx), (3.10) 

Now, consider the case of finite I and r. In this case, (3.8), (3.9) and (3.10) define a regular 
Sturm-Liouville problem. It is well-known (see for example [33,34]) that the eigenvalues A„ of 
all regular Sturm-Liouville problems are discrete and the eigenfunctions ^n (x) form a complete 
orthonormal basis in L^{Ij s). 

The situation is somewhat more complicated if cither I, r or both are infinite. In this case (3.8), 
(3.9) and (3.10) define a singular Sturm-Liouville problem. In general, the spectrum of a singular 
Sturm-Liouville problem may be discrete, continuous or mixed. Additionally, the eigenfunctions 
of a singular Sturm-Liouville problem may be "improper" , in the sense that they do not belong 
to L^{I, s). Nevertheless, the eigenfunctions may still be used as a complete set of basis functions 
in the same sense that je**^^ : k € R} can be used as basis functions of a Fourier transform (see 
for example p. 161 of [31] or p. 318 or [18]). 

For a general second order linear operator C = a{x)d'^^ + b{x)dx + c{x) on some interval / 
(possibly finite, infinite or semi-infinite), there exist sufficient conditions that one may check in 
order to classify the spectrum of the operator (see Chapter 22 of [34]). For ^£(°^), the linear 
second order operator considered in Theorem 3.1, a direct computation reveals that the spectrum 
is continuous when considered on infinite or semi-infinite intervals and the eigenfunctions - while 
improper - are complete in the appropriate Hilbert space. 

In the following Proposition we present an explicit solution -^ ^Pq (x), Aq > to the above Sturm- 
Liouville problem. 



Proposition 3.2. Depending on the interval I = (I, r), the solution to the Sturm- Liouville problem 
defined by (3.8), (3.9) and (3.10) — or equivalently (3.3) — is as follows: 



(a) — CO < I < r < oo, 



(6) — oo = / < r = oo, 



(c) — oo = I < r < 



(d) — oo < I < r = oo, 



where c is defined in Theorem 3.1. 



vl/(°)(x)=e— 
2 



-? 



sin {an{x — I)) 



a- 



(c' + a' ) , 



■ exp {ivx) 



A(") = -^(c2 + .^ 



(0)('-r.^ - p-=^ 



(c^ + , 



A(o) 
v&W(x) 



Ai°) = -y(c^ + ^ 



j/(a; — r)) 



j^(a; — I)) 



tn — 

r — 


7Ao.iiJ 


n e N, 


(3.12) 




(3.13) 


i^eR, 


(3.14) 




(3.15) 


V e R+, 


(3.16) 



1^ e 



Proof A direct calculation shows that the above eigenvalues and eigenfunctions satisfy (3.8), (3.9) 
and (3.10). One can easily verify that the eigenfunctions form a complete orthonormal basis in 
the corresponding Hilbert spaces. D 

Having found explicit expressions for vp^ '(^x) and Xq , and having established the completeness 
of < ^q (x) > in L^ (/, s) we are now able to present our solution to (3.7), (3.8) and (3.9). 

Theorem 3.3. Suppose 

(^(^) , ^(1) *10)) ^ - C(i) (m, n) I{„^„} + i?(i) (n) ,5„,„ , (discrete spectrum) 

U(o)^j^ii) ^^{o)\ ^ C(i)(u;,iy)l{^-^^} +£>(!' (i/) (5 (cj - j/). (continuous spectrum) 

Then the solution to equation (3.7) with BC's (3.8) and (3.9) is given by 



m 

A«=i?W(n) 
^(^Hx) ^ I al]l¥S\x)du;, 
A«=i?«(-) 



".™ (0) (0) ^{m=^n}, 

(discrete spectrum) (3-17) 

"'i^,^ ^ ,(0) ,(0) ^{"7^1^}' 

(continuous spectrum) (3.18) 



Proof. By Theorem 3.1 the spectrum of (/I*^"-') is either discrete or absolutely continuous and the 
eigenfunctions form a complete orthonormal basis in L^(I, s). We consider the discrete spectrum 
case. For every n G N, the function \l/„ (x) may be expressed as a linear combination of basis 
functions 



<^^-E<in"'- 



(3.19) 



Inserting (3.19) into (3.7) yields 



Multiplying both sides by ^m (x)s(x) and integrating with respeet to x we find 

(^(0)^^(1) ^(0))^ ^ ;^ 41) (,(0) _ ,(0)) (^(0)^ ^(0))^ ^ ,^^^ 
* A; 

C^^Hm,n)I{^^^^+D^^\n)S^.n = 4^^ (aI^^ - A^') + AW,5™,„. (3.20) 

Equation (3.20) is satisfied for all m and n by choosing a„,rn and A„ as in (3.17). Note that 
BC's (3.8) and (3.9) are satisfied by construction. The proof in the continuous spectrum case is 
analogous. D 

3.2 Option Prices 

We have now found expressions for the approximate eigenvalues A^ ss Aq + \/eXq and approx- 
imate eigenfunctions \I'^(a;,2/) « "^q (x) + ^/e^g (x) of eigenvalue problem (2.11), (2.12) and 
(2.13). We now use these expressions to specify the approximate price u'^{t,x,y) « u^^\t,x) + 
y/€U^^'{t, x) of an option. This serves as the main result of our work. 

Theorem 3.4. The approximate price of an option is given by 

u'{t,x,y) « u^°\t,x) + ^/eu'^^\t,x), 
where 



X^n ^" 5" ^" (discrete spectrum) 
/ AI 'gl '^/l 'dv (continuous spectrum) 



" - \ r .(0) (O).-r.(O) , / ,. , X ' I'^-^-'-j 

and 



E„ ( Ai^^g^^^^ + Al°'5i')*l°' + Ai°)gl°)*li) ) (discrete spect 



truni 



''^'' ^ J (AV^gl°'*l°' + A(°\9l^)*(0) + Al°'.gl")*l^)) 1 (continuous spectrum) ' ^^'^^^ 



Here, 



and 



ff (°) (i) = exp (a(°) t) , g(i) (t) = (a W i) exp (a ^ t) , (3.23) 



4°'-(*f'/^)^' (3.24) 

4''=-E^" (*"^*™) ' (discrete spectrum) (3.25) 



AW =- /"a£") f*l"\*i,^)') . (continuous spectrum) (3.26) 

The C (e") eigenfunctions '^q (x) and eigenvalues Xq are given in Proposition 3.2, and their 
O (s^') corrections ^q (x) and Xq are given in Theorem 3.3. 
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Proof. Consider the spectral representation of u'^{t, x, y) given by (2.14). We expand A'i and gg{t) 
in powers of -y/e 

Inserting these expansions as weh as the expansion for ^^ into (2.14) and collecting terms of like- 
powers of Ve yields (3.21) at O (e°)and (3.22) at O (e^/^). The expressions (3.23) are obtained 
from (2.14) by performing a Taylor series of (7^(i) about -y/e — 0. Expressions in (3.24), (3.25) and 
(3.26) can be obtained from the BC u'^{0,x,y) — h{x). We make the choice u(°)(0,x) = h{x) and 
u^^\Q^x) — 0, which is consistent with the choice made in [10]. Temporarily specializing to the 
discrete spectrum case we note 

.(°)(0,x) ^ h{x) ^ 5:4°)*(0)(,) ^ ($(0),;,) ^ ^^(0) (*(0), vpW)^ = 40), 

m m 

Likewise 



n(i)(0,x) = = 5: (A«*(;])(-) + At^^^^Hx) 

m 



Ti ' / J rn \ n ^ rn 






The continuous spectrum case is analogous. D 

Corollary 3.5. The function ^.u^^'{t^x) is linear in the group parameters 

V,^:^V'^^{m^V^V2, Vi:^-Ve^{f<l^')^V^V3. (3.27) 

Proof. From (3.6) we see that the operator A^''^' is linear in V2 and V3. By Theorem 3.3 it is clear 
that ^q (x), a,p,q, and \q are linear in V2 and V3 as are gq (t) and Aq by (3.23), (3.25) and 
(3.26). Finally, because V2 and V3 do not appear in *g (x), Xq and gq it is clear from (3.22) 
that u'^^\t,x) is linear in V2 and V3. Thus, ^/eu''^\t,x) is linear in V2 and Kf. D 

3.3 Equivalence to Black-Scholes and to Fouque-Papanicolaou-Sircar 

[10] 

In this section, we will show that u'"^ (t, x) corresponds to the Black-Scholes price of an option 
with Black-Scholes volatility equal to vct^. We will also show that u'^^^(t, x), the 0{^/e) correction 
to u^'^' (t, x) due to fast mean-reversion of the volatility, is the same correction as that obtained 
in [10]. This equivalence relation will enable us to establish the accuracy of the pricing approxi- 
mation u'^{t,x,y) ss u^^'{t,x) + y/eu^^'{t,x) for the case of European options. 

Theorem 3.6. Let u^^{t,x) be the Black-Scholes price of an option with with payoff (2.1) and 
let the underlying have volatility \/a^. Then 



u^S 



{t,x) =u(0)(i,x). 
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Proof. In the Black-Scholes model, the underlying is assumed to follow geometric Brownian motion 
with risk- neutral drift /i and volatility V(t^. The Black-Scholes price u^^ {t, x) of a an option with 
payoff (2.1) solves the following PDE with BC's 

0^i-dt+Cx)u^'^, {t,x)e[0,T]xI, (3.28) 

h{x)^u^^{0,x), (3.29) 

O^u^'^itJ) if />-oo, (3.30) 

= u^'='{t,r) if r<oo, (3.31) 

where 



Cx^i^^i--a'jd, + -cT'd',,. (3.32) 

By construction u(°)(t,x) satisfies BC's (3.29), (3.30) and (3.31). Hence, by the uniqueness of 
the solution to the above linear PDE problem, in order to establish the equivalence of u'"-* {t, x) 
to u^^{t,x) we need to show that u("^(i,x) satisfies PDE (3.28). To this end we note that 
Cx = (£'-"■'). Now, specializing to the discrete spectrum case, we see that 

n n 



The calculation in continuous spectrum case is analogous. Hence, we deduce that u {t,x) — 
u^°\t,x). D 

Theorem 3.6 is consistent with the findings of [10], where it was found that the 0{e'^) price of 
an option was given exactly by u^^{t,x). 

Theorem 3.7. Let y/eu^^^{t,x) be the O (e^^^) correction to the Black-Scholes price u^^{t,x) 
of an option with payoff (2.1) as calculated in [10]. Then 

u^P%t,x)=u^^Ht,x). 

Proof. It is established in [10] that the EPS correction u^^^{t, x) to the Black-Scholes price of an 
option u^^{t,x) satisfies the following PDE and BC's 

= «^^^(0,x), 
= u^P^it,l) 
O^u'^P'it^r) 

By construction u'-^^t^x) satisfies BC's (3.34), (3.35) and (3.36). Hence, by the uniqueness of the 
solution to the above linear PDE, in order to establish the equivalence of u(^'(t,x) to u^^^{t,x) 
we need to show that u'^^^t.x) satisfies PDE (3.33). Using u^^{t,x) = M(")(t,x), Cx = {C^"'>), 
expression (3.21) for u'''^\t,x), expression (3.22) for u'-^\t,x), a straightforward but tedious cal- 
culations yields (in the discrete spectrum case) 



(t,x)e[o,r]x/. 


(3.33) 




(3.34) 


if / > — oo. 


(3.35) 


if r < oo. 


(3.36) 



a..(^)=E4°)^)AW(l+tA("))vI;(o), 

n 
n 
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Inserting the above equations into (3.33) verifies that u'-^\t,x) satisfies PDE (3.33). The calcu- 
lation in the continuous spectrum case is analogous. Hence, we deduce u^^'{t,x) ~ u^^^ {t,x). 

D 

Conveniently, the equivalence relation 

M(°)(i,x) + Viu(i)(i,x)=M^^(t,a;) + Viw^^^(i,x), 

establishes the accuracy of our pricing approximation. 

Theorem 3.8. Under assumptions 1 - 4 of section 2.1 and under the assumption of bounded f{y) 
we have the following accuracy results: 

1. For European options with smooth and hounded payoffs, for all t < oo and for x, y € R 

u'{t,x,y) - (u'^^'Ht.x) + V^w(i)(i,x)) I = 0{e). 

2. For European call options, for all t < oo and for a;, y € M 

u^(t,x,y)-(M(o)(t,x) + Viu(i)(i,x)) =0(e|loge|). 



Proof. The proofs of 1 and 2 are given in [10] and [12] respectively. The proof for unbounded f{y) 
satisfying condition 5 of section 2.1 can be found in [13]. D 

We remark that the accuracy results of Theorem 3.8 are valid when e is much smaller than the 
life of the option. The reason for this is that our pricing approximation depends on the process 
y/ having sufficient time for the time-average of f(Y^) to approach its ensemble average 

i J^ f (r/) ds^ij^ f (n%) ds^l- f^'^ f (yj) du ^ if') . (3.37) 

The accuracy results of Theorem 3.8 are for fixed t. It is clear from (3.37) that convergence is not 
uniform in t. For barrier options, if x is near an cndpoint I or r, the life of the option may be 
of order e due to X hitting a barrier prior to the time of maturity T. Thus, convergence is not 
uniform in x. A detailed analysis of the accuracy of our pricing approximation when x is near 
an endpoint would require boundary layer analysis. Such an analysis is beyond the scope of this 
paper. 

4 Practical Implementation 

In this section we discuss the practical implementation of our methods. In sections 4.1, 4.2 and 
4.3 we provide three examples, which show how the results of Sections 3.1 and 3.2 can be used 
to specify the price of an option. In section 4.4 we sketch how our results can be extended to 
price rebate and knock-in options. And, in section 4.5 we provide a recipe for calibrating the fast 
mean-reverting class of models to the market using European call option data. 

4.1 Example: European Call Option 

The payoff of a European call option with strike price K — e and time to maturity t can be 
expressed in the framework of (2.1) by choosing ^ 

h{x) ^ {e'' - e'')^ , /= (-00,00). 



•^Notc that the payoflf h{x) is not in L^{I, s). This can be dealt with by appealing to the theory of generalized 
Fourier transforms. We discuss this further when we calculate A), ' . 
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Note that as X can not leave / — (—00, 00) in finite time we have t = 00 and !{,-> (} = 1. Hence, 
the payoff of the option is given simply by h (Xt), which is as it should be for a European option. 
To calculate the approximate price of a European call option u''^\t^x) + y/eu'^^\t^x) the first 
thing we must do is find expressions for the approximate eigenfunctions ^1/ (x) + ^'^l, {x) 
and eigenvalues \\^ + \/e\y ' . The O (e°) eigenfunctions W^ {x) and eigenvalues Xi, are given 
explicitly by (3.13) and (3.14) of Proposition 3.2. To find the O (e^^^^) corrections vJ/J^ ' (x) and 
Ai/ we use Theorem 3.3. We note 

(vl/W,^W vI/W)^ = CW(c.,i.)I{^^,j +i?W(;.)(5(c. - ^.), 

D'^^\v) - VsP, + V2C,, Cu = (iv - cf - (iv - c). 

Hence, from (3.18) we find 

*l''(a;) - 0, aW = 0, A(,i) = V^fi, + V^C,,. 

We must now find expressions for gl, (t), g\, (t), J\i and J\i, . This can be accomplished using 
Theorem 3.4. Having identified Xi and A^, , we read gi, (t) and gl (t) directly from (3.23). The 
coefficients A^ and AI are obtained from (3.24) and (3.26). We have 

— exp (ivx) (e^ - e'') -^e^^'-^dx (4.1) 

47r ^ / 0-2 




/fj27r {iv - c- l){ii'- c)' 

Note that integral (4.1) will not converge unless we impose Im [v\ > (c+1). Thus, in deriving result 
(4.2), we have implicitly assumed u = Vr -\- ivi and fixed vi > {c+ 1). The process of extending 
the domain of the variable of integration into the complex plane, which is contained in the theory 
of of generalized Fourier transforms [32], enables us to extend our results to European options 
with payoffs h{x) ^ L'^{I,s). It is important to note, however, that because of the condition 
Vi > (c+ 1), when evaluating integrals (4.3) and (4.4) below we must make sure to set v — Vr + ivi 
and dv ~ dvr (i.e. integrate over a contour parallel to the real axis in the complex plane). 

Having found Ai, and Ai, , the approximate option price u'^{t,x,y) « u^'^'(t,x) + y/eu^^'{t,x) 
can be found from (3.21) and (3.22). We have 

/oo 
A^^'>g^J'Ht)¥^Hx)dv, (4.3) 

-00 
/OO 
A^^'>gi'Ht)¥^Hx)dv. (4.4) 

-00 

Now, recall from Theorem 3.6 that u*^°^(i,a;) — u^^{t,x), the Black-Scholes price of a Euro- 
pean option with volatility vct^. And recall from Theorem 3.7 that u'^^\t^x) = u^^^{t,x), the 
correction to the Black-Scholes price due to fast mean-reversion of the volatility, as calculated 
in [10]. Finally, recall from Corollary 3.5 and that u^^'{t,x) is linear in the group parameters 
parameters V2 and V^, defined in (3.27). For European call options, it was shown in [10] that 
the group parameters parameters V2 and V^ have a very specific affect on the implied volatility 
surface induced by fast mean-reverting stochastic volatility models; a change in V2 corresponds 
to an adjustment of the overall level of implied volatility and a change in V^ corresponds to an 
adjustment of the at-the-money skew. This structure leads to a remarkably simple calibration 
procedure, which we outline in section 4.5. The effect of V2 and Kf on European call prices 
and the corresponding effect on the implied volatility surface is demonstrated in figures 1 and 2 
respectively. 
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Figure 1: Prices of European call options are plotted as a function of St, the current price of the 
underlying. In these sub-figures, t = 1/2, /i = 0.05, vo^ = 0.34 and exp(fc) = 2. In sub-figure (a), 
we set V^ — 0, and vary V2 from —0.03 (red, dot-dashed) to 0.03 (blue, dashed). In sub-figure 
(b), we set V2 = 0, and vary V^ from —0.01 (red, dot-dashed) to 0.01 (blue, dashed). In both 
sub- figures the solid line corresponds to the Black-Scholes price of the option (i.e. V2 =V^ = 0). 
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(a) 



(b) 



Figure 2: Implied volatilities of European call options are plotted as a function of strike price K. 
In these sub-figures, i = 1/2, /i = 0.05, Vct^ = 0.34 and St — 2. In sub-figure (a), we set V^ — 0, 
and vary V2 from —0.01 (red, dot-dashed) to 0.01 (blue, dashed). In sub-figure (b), we set V2 = 0, 
and vary V^ from —0.01 (red, dot-dashed) to 0.01 (blue, dashed). In both sub-figures the solid 
line corresponds to / = ct (i.e. V2 — V^ = 0). 

4.2 Example: Up-and-Out Call Option 

The payoff of an up-and-out call option with knock-out barrier R = e'^ < 00, strike price K = 
e^ < e^ and time to maturity t can be expressed in the framework of (2.1) by choosing 



h{x) 



Hxei}i 



I = (— oo,r) . 



To calculate the approximate price u(°^(i, a;) -I- y/eu'^'^\t,x) of such an option we must first find 



expressions for the approximate eigenfunctions "^t'' [x) + \/e "^V^' {x) and eigenvalues \\^'' + ^/e Xl . 



,(0), 



,(!)/ 



>(0) 



The O (e°) eigenfunctions "$i (x) and eigenvalues AJ^ are given explicitly by (3.15) and (3.16) 
of Proposition 3.2. The O (e^^^) corrections ^Jy (x) and Xl are found using Theorem 3.3. We 
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calculate 

X = 2c + 1, r]^=.u^ - (3c2 + 2c) , 

e. = -t^' + (c' + c) , 7, = (3c + ly - (c3 + c2) . 

Now, from (3.18) we have 

^i'\^) = / 4!i.*i"'(^)rf'^, aW = -— -3 (y^X + nv.) , (4.5) 

Jo cr TT (cj^ — l^-^j 

aW = V2C. + V3J.. 

To find gl (i), gj) (i), Al and Ajy we use Theorem 3.4. Having identified Xl and Xl , the 
functions gl (i) and gi/ (i) are read directly from (3.23). The coefficients Ai ' and Ai ' are 
obtained from (3.24) and (3.26) respectively. We have 



AW^J^^e 



c2 + i/2 (1 + c)2 + 1/2 



4 fi/e^+'^'' (xcos(j/(r - A:)) + ^j, sin (j/ (r - fc))) 



cr27r V (c2 + l/2)((l + c)2 + t/2) 








Finally, the approximate option price u'^{t, x, y) « u^^'{t, x) + y/eu^^'{t, x) can be found from (3.21) 
and (3.22). 



00 /-OO 



u^'\t,x)= A(0)g(0)(i)vl/(0)(x)d^, 

w(i)(t,x)= / AWgW(i)vl/W(x)d;. 

g(^\t) (A(°)aW *(")(x) - 4") ««*(")(,)) d^rf.. (4.6) 

Note that, while the double integral in (4.6) is finite, it blows up along the line uj = 1^. This 
complicates numerical integration schemes. A method of dealing with this issue is provided in 
appendix A. Figure 3 demonstrates the effect of parameters V2 and V^ on the price of an up-and- 
out call option. 

4.3 Example: Double-Barrier Knock-Out Call Option 

The payoff of a double-barrier knock-out call option with barriers L — e^ and R = e^, strike price 
K = e^ (with —00 < I < k < r < 00) and time to maturity t can be expressed in the framework 
of (2.1) by choosing 
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(a) 




(b) 



Figure 3: Prices of up-and-out call options are plotted as a function of St, the current price of the 
underlying. In these sub-figures, t = 1/12, /i = 0.05, vo^ = 0.34, exp(fc) = 2 and exp(r) = 2.5. In 
sub-figure (a), we set V^ = 0, and vary V2 from —0.01 (red, dot-dashed) to 0.01 (blue, dashed). 
In sub-figure (b), we set V2 = 0, and vary V^ from —0.001 (red, dot-dashed) to 0.001 (blue, 
dashed). In both sub-figures the solid line corresponds to the Black-Scholes price of the option 
(i.e. V^^ = V,^ = 0). 



To calculate the approximate price u^^'{t,x) + \/eu^^'{t,x) of such an option we must first find 
expressions for the approximate eigenfunctions "iii ' (x) + ^/e"^i ' {x) and eigenvalues Xl' + ^/eXl'. 
The O (e°) eigenfunctions ^„ (x) and eigenvalues An are given explicitly by (3.11) and (3.12). 
The O (e^^^) corrections ^„ '{x) and A„ are found using Theorem 3.3. We calculate 



^vI.W,^(i)vI,Wj^ = C(i)(m,n)I{„,^„}+i?(i)(n)5,„,„, 

X = 2c+1, 
7j„ = al - (3c2 + 2c) , 

£.n = -al + (c^ + c) , 



(T^2C„ + V37J„) , 



7„-(3c+l)«2-(c3 + c2). 



Now, from (3.17) we find 






.(0)/ 



,(1)/ 






l(0) 



^ ((-1)™+" - 1) 



Z(yiji0^i2 



[r - /) {al 



{V2i + 1^3^n) , 



i(i) 



(0) 



i(l) 



In order to find gn (t), gii (t), An and An ' we use Theorem 3.4. Having identified A„ and A 
gn (t) and gn (t) are read directly from (3.23). The coefficients An and An are obtained from 
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(3.24) and (3.25). We have 



Al°) = 



4 / e" 

-(-l)"e^'-a. 



4 ie''+''''{xanCOs{an{k~l)) - f„ sin(a„(fc - /))) 



a2(r-0 



(c2+a2)((l + c)2 + a2) 



/ J in ni,n' 
m 

Finally, the approximate option price u'^{t, x, y) « u^^'{t, x) + y/eu^^' (t, x) can be found from (3.21) 
and (3.22). 

B(°)(t,x)=E4"^gi°)(t)*(")(.), 

n 

z.W(t,x) = E4"^ffW(0vi'(")(x) 

n 

n m 

Figure 4 demonstrates the effect of the parameters V2 and V^ on the price of a double-barrier call 
option. 
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Figure 4: Prices of double-barrier knock-out call options are plotted as a function of S't, the 
current price of the underlying. In these sub-figures, t = 1/12, /i = 0.05, \fa^ = 0.34, exp(fc) = 2, 
exp(/) — 1.5 and exp(r) — 2.5. In sub-figure (a), we set V^ — 0, and vary V2 from —0.01 (red, 
dot-dashed) to 0.01 (blue, dashed). In sub-figure (b), we set V2 = 0, and vary V^ from —0.001 
(red, dot-dashed) to 0.001 (blue, dashed). In both sub-figures the solid line corresponds to the 
Black-Scholes price of the option (i.e. V2 = V^ = 0). 



4.4 Brief Note on Knock-in and Rebate Options 

To this point, we have considered only options with payoffs given by (2.1). Options that fit within 
this framework include European options and knock-out style options. In fact, our pricing results 
can be extended to include knock-in and rebate style options as well. The focus of this section 
is to give an idea of how this extension can be done. For the sake of brevity, the proofs in this 
section will be kept short and will contain only the main ideas needed for the full proofs. 
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First, we consider a knock-in option. Such an option has a payoff of the form 

PayoSknock-tn = H^t) I{r<T}, T = inf{t > : Xj <^ /}, 

/ := (Z,r), — oo <l<r< oo, 

We compare this to a knock-out option and a European option (both of which we have already 
priced) whose payoffs can be reformulated as follows 

Payoff;,„„^;,_„„t = hiXr) I{r>T} Pl^yof^European = H^t) ■ 

We note that the knock-out option has payoff h{XT) on the event {t > T} whereas the knock-in 
option has payoff h{XT) on the event {r < T}. The European option has payoff h{XT) regardless 
of when r occurs. We can use this information to specify the price of a knock-in option. Briefly, 

Payofffe„ocfe-o«t + Payofffe„ocfe-m = K^t) I{r>T} + KXt) I{r< T} = KXt) = ^^yO^European- 

Taking expectations on both sides, it follows that the price of a knock-in option is just the price 
of a European option minus the price of a knock-out option. A more detailed discussion of the 
knock- in knock-out parity relation can be found in [4]. 

Now, consider a double-barrier rebate option (the the single-barrier case is analogous). The 
payoff of a double-barrier rebate option is given by 

Payoff,,fc,te ^h{Xr), t = inf{i > : Xj ^ /} A T, 

/ := {I, r), — oo < I < r < oo, 

h:lU{l}U{r}^R+, h(l) ^ Ri, h{r) ^ Rr. 

Note that the restriction h(l) = h(r) = of equation (2.1) has been relaxed. The payoff of the 
above option can be described as follows: if the log of the underlying does not exit {I, r) prior to 
time T, the option has payoff h{XT), otherwise the option pays a rebate i?; > if the log of the 
underlying exits at I or pays a rebate Rr > if the log of the underlying exits at r. 

The price P^ of such an option at time s < T is given by (2.2). However, the first term in 
(2.2) is no longer zero, due to a rebate being paid at time r on the set {r < s}. Additionally, the 
function P'^{s,x,y) now satisfies the following PDE and BC's (see Chapter 9 of [29]) 

0={^s-^l + H'x.y) Prebate, («, X, y) G [0, T] X / X K, 

h{x) ^ P^^^^t^{T,x,y), 

Rl = Prebatei^i '7 2/); 

where we have added the subscript rebate to indicate that we are specifically considering rebate 
options. In terms of u^^fj^^^{t,x,y), whose relation to ^rehate (*'■*"' 2/) i^ defined in (2.6), we have 

{t,x,y)e[0,T]xIxR, (4.7) 

(4.8) 

(4.9) 

(4.10) 

In order to specify the approximate price (^u'-^' + y/eu^^^) . of a rebate option, we shall need 
the following Lemma. 

Lemma 4.1. The price u^ebatei^T^^y) '^f '^ rebate option can be expressed as 

Kebateit, X, v) = e^''^^\x, y) + v'{t, X, y), (4.11) 
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= 


-^{-dt+C'x.yX 


'Mate'' 


h{x)^ 


= debate (0,2 


■.y), 




ef^'Ri = 


= ^rebate\tjh 


y). 




e^''Rr = 


- "rebate (*,''' 


,y)- 





where ^'^{x,y) satisfies 

0-(-M + /:5f,y)*', ix,y)eIxR, 

i?, = r(r,y), 

and v'^{t,x,y) satisfies 

= (-a* + r^f^y) w% (i,x,y)e [0,cx))x/xR, (4.12) 

/i(x)-$^(x,y)=«^(0,x,y), (4.13) 

= «^(t,/,y), (4.14) 

0^v'-{t,r,y). (4.15) 

Proo/. The proof is by substituting (4.11) into (4.7), (4.8), (4.9) and (4.10). D 

Note the similarity of equations (4.12), (4.13), (4.14) and (4.15) to equations (2.7), (2.8), (2.9) 
and (2.10); the only difference is that f '(i, x, y) has an e-dcpendent BC in (4.13) whereas u'^{t, x, y) 
in (2.8) does not. The e-dependent BC requires a minor modification of the asymptotic analysis 
of section 3. The result of this modification is contained in the following Theorem. 

Theorem 4.2. The approximate price {^-^'{t^x) + ^u^^'{t,x)^ of a rebate option is given 

by 

^i(°) (t, x) + Viu(i) (i, x)\ = e^* ($(") (x) + V^ $(1) (x)) + (i;(") (t, x) + ^fev^^^^ (t, x)\ 

where ^^^'{x) and ^^^'{x) satisfy 

0= (/i-(£(o) ))$(«), i?, = <i>(0)(0, i?. = $(°)(r), 

^(i)$(o) = r^_/£(o)\U(i), = $«(0, = $(i)(r). 



T/ie functions v^^> {t, x) and v^^' (i, x) have spectral expansions given by the right hand side o/ (3.21) 
and (3.22) in Theorem 3.4 where An and An are now given by 

41) =_(.,(0)^ <,(!)) -^A(^)(v,(0),V,(l) 



Proof The proof follows from Lemma 4.1 and by modifying the analysis of section 3 to account 
for the e-dependent BC from (4.13). D 



We note that <i>'^°^(x) is given by 



(0) /i?^e-sinl^-0(i + ^))+i?.e^'sinh((r-x)(i + i^)) 



sinh((r-0(i + ^)) ; 

We omit the expression for $'^^'(x) for the sake of brevity. 

4.5 Calibration 

In this section we will briefiy discuss how one can calibrate the class of fast mean-reverting models 
to the market using European call option data. 

One of the great advantages of the option-pricing framework developed in this paper is that, 
although the fast mean-reverting volatility process adds five parameters {y, e, v, p, y) and two 
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unspecifiod functions (/ and A) to the Black-Scholes framework, specific knowledge of these pa- 
rameters and functions is not needed to specify the approximate price of an option. Instead, the 
parameters and functions listed above are replaced two group parameters, V2 and V^', given by 
(3.27). What is more, V2 and V^ are defined consistently throughout this paper irrespective of 
the type of options being considered. That is, the group parameters V2 and V^^ that are used to 
give the approximate price of a European call option are the same parameters that are used to 
give the approximate price of e.g. a double-barrier knock-out option. Thus, one can use (liquid) 
European call option data to calibrate the class of fast mean-reverting stochastic volatility models 
to the market. Once this is done, the obtained group parameters can be used to price (illiquid) 
exotic options. The following calibration procedure is suggested in [10]: 

1. Using (liquid) European call options, fit observed implied volatilities /^ as an affine function 
of the log moneyness to maturity ratio (LMMR^j ) 

/,, = b + aLMM%, LMM% = \ogiK,,/St)/iT, - t), 

where lij is defined implicitly through 

u^^ (r„ K,j , /,, ) = u^'""^"* (T, , K,, ) . 

2. The group parameters V2 and V^^ are then given by solving 



We note that a, the average level of volatility of the underlying, which can be obtained from 
historical returns data, is needed to determine ¥2- 

3. Use the obtained values for a, V2 and V^ to give approximate prices for (illiquid) exotic 
options. 

The above calibration scheme was tested with single-barrier knock-out options in the context of 
credit risk in [14], where it was shown to work well. 

5 Conclusion 

Using elements from spectral analysis and singular perturbation theory, we have presented a 
systematic way to obtain the approximate price of a variety of European and path-dependent 
options in a fast mean-reverting stochastic volatility setting. One key feature of our technique is 
that we were able to maintain correlation between the stock-price and volatility processes via two 
Brownian motions and still produce pricing formulas for double-barrier options. To our knowledge, 
this is the first paper to address this issue. Extending our techniques to more sophisticated models 
is an on-going process. A logical next step, for example, would be to add a fast mean-reverting 
factor of volatility to a model such as CEV or Heston as done in [9] or to add a slow- varying factor 
of volatility to the class of models considered in this paper. 
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A Addressing Numerical Integration Difficulties 

In this section we demonstrate how to accurately evaluate the double integral in equation (4. 
which we repeat here for clarity 



OO /"OO 



J- / 3(°)(4")aW*(")-AL0)aW*(°))rf^d.. (A.l) 



'0 "'0 

The difficulty in numerically evaluating (A.l) is that, for most Al , the integrand blows up as 
w — )• J^. This is due to the factor of 1/ (i^^ — w^) which appears in al.l, (refer to equation (4.5) for 
details). Thus, as it is written in equation (A.l), numerically evaluating J would require adding 
and subtracting some very large numbers, which most numerical integrators are not very well- 
equipped to do. Thankfully, there are a few numerical tricks we can perform in order to facilitate 
numerical evaluation of (A.l). To begin, we establish some notation. Let 

and make the following change of variables 

'^iu,v) ==—=(«-«), 
w(w,w) = — =(m + w). 



Now, we define 



so that 



G{u,v) -.^ giv{u,v)), 
H{u,v) :~ h{i'{u,v),uj{u,v)), 



J = I I G{u,v)H{u,v)dvdu. (A.2) 



So far, everything we have done is cosmetic; the integrand of equation (A.2) still blows up near 
V = (which corresponds to i/ = lu). Note, however, that H{u, v) = —H{u, —v). As such, we may 
write equation (A.2) as 

J= / Hiu,v)iG{u,v)-Giu,-v))dvdu. (A.3) 

Jo Jo 

The integrand in equation (A.3) is well-behaved throughout its domain. Figure 5 illustrates how 
this simple trick smooths out the singularity. 
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